Time resolved fission in metal clusters 



O 

o 

(N 

> 
O 



5h 
43 



o3 



o 
o 



> ■ 

00 ' 



-xt- : 
o . 

-i— > ■ 
a ■ 

E 

T3 ■ 

c : 
o . 
o ■ 



X 



P. M. Dinh 1 , P.-G. Reinhard 2 , E. Suraud 1 , 2 
1 Laboratoire de Physique Theorique, Universite Paul Sabatier, 
118 Route de Narbonne, F-31062 Toulouse, Cedex 04, France and 
2 Institut fur Theoretische Physik, Universitat Erlangen Staudtstr. 7, D-91058 Erlangen, Germany 

We explore from a theoretical point of view pump and probe (P&P) analysis for fission of metal 
cluster where probe pulses are generalized to allow for scanning various frequencies. We show that 
it is possible to measure the time the system needs to develop to scission. This is achieved by a 
proper choice of both delay and frequency of the probe pulse. A more detailed analysis even allows 
to access the various intermediate stages of the fission process. 



Pump and probe (P&P) experiments have become a 
widely used, versatile and powerful tool for time resolved 
studies of dynamical reactions in all sorts of electronic 
physics from molecules 0, Q to bulk [j| . Experimental 
and theoretical investigations of P&P analysis for clus- 
ters as intermediate size systems have started more re- 
cently and are mostly still on an exploratory level 
see also [?| for an overview. Clusters are here a great 
challenge because they are much more complex than 
molecules and there is thus a much larger variety of pos- 
sible scenarios. This calls for theoretical explorations and 
there exist already several of them, see Q. In the spirit 
of the studies in 0,01, we propose to exploit the Mie plas- 
mon resonance as a handle to track the time evolution of 
the cluster shape. A most dramatic evolution of shapes 
is experienced in cluster fission 

[HIIll, see also for 
a review. The various shapes along the fission path are 

related to characteristic resonance spectra which pro- 
vides a unique opportunity for P&P analysis of the as- 
sociated time scales. As we shall see below a standard 
P&P scenario, augmented by a frequency scan, indeed 
allows to access a fully time resolved analysis of fission 
dynamics. This is an interesting task as such and it car- 
ries farther reaching information. The time the system 
needs to attain scission (namely to actually split into two 
pieces) provides an indirect but crucial clue to the vis- 
cosity of clusters . Fission is also a well known prop- 
erty in another class of finite fermion systems, namely 
atomic nuclei. It is actually by studying fission dynamics 
in atomic nuclei that the first experimental estimate of 
nuclear viscosity could be attained (ljj . 

The dipole plasmon plays a key role in the dynamics 
of metal clusters, both in the linear and non-linear do- 
mains 0, flil . It will again provide here the handle for 
the P&P analysis. The theoretical modeling needs to ac- 
count properly for both, electronic and ionic degrees of 
freedom. This is achieved by the Time Dependent Local 
Density Approximation (TDLDA) for electrons, coupled 
to Molecular Dynamics (MD) for ions, the coupling be- 
ing achieved by a local pseudo-potential 01- Absorbing 
boundary conditions are used throughout, which allows 
a proper description of ionization. We use a cylindri- 
cally averaged description of the electrons to simplify the 



extensive calculations Details on the TDLDA-MD 
approach can be found in 0, 0] ■ 

The laser field is described in the dipole approxi- 
mation as a time-dependent external potential Vj as — 
Eo-d/i as (t) cos(wt), with the dipole operator d, the laser 
polarization and amplitude Eq, and fu s (t) the temporal 
profile of the laser pulse, here chosen as a sin 2 profile. 
Pulses are short (FWHM typically less than 50 fs) as 
required for a P&P scenario. Laser intensities are kept 
moderate and tuned to provide the proper pump excita- 
tion and the proper probe analysis in terms of ionization. 
Pump as well as probe pulses have the same polarization 
but different frequencies and FWHM (see below). 

Both electronic and ionic observables are recorded as 
a function of time. The two key observables for the P&P 
scenario are ionization and the electronic dipole moment 
with respect to the ionic center of mass, which provides 
the plasmon response by Fourier transformation into the 
frequency domain, |18l Il9| . The ionization is computed 
as a function of time, by recording the number of elec- 
trons lost through the absorbing boundaries (N esc ). Ions 
are treated as classical particles by standard MD propa- 
gation of positions and momenta. The P&P analysis will 
thus provide a way, with optical techniques, to access this 
ionic dynamics experimentally. 

We consider the cluster Nai4 initially excited by a 
pump laser to a 3 + charge state. The laser is linearly 
polarized along the symmetry axis of the cluster. The 
Nai4 +++ is unstable with respect to fission an d thus 
starts to elongate until it breaks up. With the pump laser 
pulse used here (intensity: 5 x 10 10 W/cm 2 , frequency: 
2.3 eV, FWHM: 36 fs), the average fission path leads to 
Nai4 +++ — ► Nag + + Na6 ++ . A word of caution is in 
place here: TDLDA-MD provides an average description 
and does not give access to all the possible fission chan- 
nels. Still, even with that limitation TDLDA-MD suffices 
for a first exploration of P&P analysis of cluster fission. 
The fission scenario which we have in mind allows a se- 
lective identification of either the fissioning system or the 
fission fragments in terms of optical response. It is indeed 
typical that the plasmon frequencies of the fragments lie 
safely above the plasmon frequency of the elongated fis- 
sioning cluster. The key mechanism which makes the 
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FIG. 1: Schematic drawing of the evolution of the plasmon 
response of a metal cluster along its fission path. Inserts in- 
dicate typical ionic structures. The P&P analysis has to be 
considered here in terms of both delay and probe frequency 
as sketched on the figure. For details, see text. 
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analysis possible is the fact that the ionization yield is 
strongly correlated to the plasmon resonance: The closer 
the frequency of the probe laser to the plasmon frequency 
the larger the ionization and vice versa [20| . The plasmon 
response, in turn, is usually concentrated in a well defined 
frequency domain and provides a direct fingerprint of the 
underlying ionic structure. Scanning the probe ionization 
thus provides a strong clue to the ionic structure and its 
time evolution. The position (and structure) of the plas- 
mon resonance changes in a characteristic manner along 
a fission path 0|: the elongation leads to a progressive 
red shift of the plasmon resonance mode along the fission 
axis. However, as soon as the system has fissioned, the 
now more compact fission products both lead to a signifi- 
cantly blue shifted response. This is sketched in figure^ 
Because of the expected strong changes of the resonance 
landscape, the probe pulses have to hit with the proper 
frequency as well as the proper delay time. We will thus 
consider a generalization of the standard P&P scenario 
to scanning probe frequencies. 

First let us start with a standard P&P analysis by fix- 
ing the probe frequency w pro b and by varying its time 
delay. For the present test case, the Mie frequencies of 
the two fission fragments are about 0.7 to 1.5 eV above 
that of the elongated Na i4 +++ which lies around 1.9 eV. 
By choosing a probe frequency at about 3 eV, we should 
thus be able to track down the scission time by scan- 
ning the dynamics in terms of time delay. Before scission 
the fragmenting system should appear more "red" while 
afterwards the system will appear more "blue". Indeed 
the observable here is the additional ionization due to 
the probe laser with respect to the mere pump ioniza- 
tion, defined as AN CSC = N csc [probe] — A^ osc [pump] . That 
quantity depends slightly on time for a given delay. We 
consider here an "asymptotic" value, recorded 500 fs af- 
ter the probe pulse. The AiV esc reflects the closeness of 
the probe frequency to that of the plasmon. For delays 



FIG. 2: Extra ionization due to the probe pulse (as compared 
to mere pump ionization) as a function of delay time for three 
typical frequencies. One clearly spots the three regimes asso- 
ciated to the three frequencies : low frequency dominant at 
short delay (full line), medium frequency dominant for inter- 
mediate delays (dotted line) and high frequency characteristic 
of post-fission and thus large delays (dashed line). See also 
figure for completeness. 

lower than the scission time, w pro be remains far away from 
the fragment plasmon and this off-resonance gives vanish- 
ing or low extra ionizations. However as we increase the 
probe delay steadily a sudden increase of AiV csc should 
occur, corresponding to the first resonance of w pro be with 
one or both fragment plasmons. Thus the first delay at 
which we observe a significant ionization after probe laser 
gives us an estimate of the scission time. This is indeed 
what we obtained for w pro be = 3.13 eV, see the full curve 
in figure|21 the obtained extra ionization, as a function of 
delay, is insignificant up to 1 ps and jumps around 1.25 
ps. 

An alternative strategy is to use a low probe frequency 
close to that of the Na i4 +++ compound system, and to 
check the disappearance of A7V CSC with increased delay 
time. The picture can be complemented by tracking 
also an intermediate frequency. These two scenarios are 
shown together with the high frequency calculations in 
figure |21 The low frequency of 1.77 eV (long dashes) at- 
taches to the elongated and still connected cluster while 
the intermediate frequency (dotted curve) attaches the 
final stages of pre-scission. As expected the three signals 
are nicely complementing each other and allow to con- 
clude on a scission time of order 1.2 ps for the given test 
case. 

The full richness of resonance spectra at a given stage 
of fission can be studied in even more detail by fixing the 
time delay and scanning a dense variety of probe laser 
frequencies. This way of analysis is illustrated in figure 
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FIG. 3: Extra emission due to the probe pulse (as compared 
to the pump ionization) as a function of probe frequency for 
a few typical delays along the fission path. Ionic structures 
at the corresponding stages during the process are shown as 
inserts. 



which demonstrates that the plasmon peak position 
(see vertical lines at the maximum of each spectrum for 
guidance) and its fragmentation pattern change dramat- 
ically on the way to scission. In turn, a probe pulse with 
fixed frequency will come into resonance only during a 
rather well defined delay window, as was seen in figure|2] 
A world of interesting detailed information is thus con- 
tained in figure when scanning the probe frequency. 
Such multi-color P&P setups are, however, an enormous 
experimental challenge. The trend of the plasmon peak 
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FIG. 4: Time evolution of various observables related to ionic 
structure along the fission path. Bottom: Plasmon frequency 
estimated from a moment analysis of the emission distribu- 
tions (see text for details). Because of the too spread response, 
we only indicates values of Ei and E% between 300 and 800 
fs. Middle: Transverse mean square radius of the cluster be- 
fore 1.4 ps and the two fission fragments for later times. Top: 
Total ionic mean square radius. 



sketched in figure ^ complies with figure For 500 fs 
delay, the spectrum is too fragmented to define unam- 
bigously a plasmon peak. This holds for the transitional 
interval 300-800 fs, as hinted by a dotted regime in fig- 
ure n 

Thus tuning probe delay time and frequency comple- 
ment each other to analyze the fission process in great 
detail. A summary analysis can be performed in terms 
of moments of the resonance spectra [2lJ. To that end, 
we compute weighted averages over the ionization spec- 
tra E n = ^jjdu; AN QSC uj n j J dw AiV 0SC uj n - 2 . Good ref- 
erence averages are E\ and E% which put moderate em- 
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phasis on the lower (Ei) and the upper (£3) side of the 
spectrum. A well peaked resonance response is signaled 
by E3 ~ E±. The difference between the two values in- 
dicates the width of the spectrum. The results for E\ 
and E3 are shown in the bottom panel of figure 0] As 
can also be seen from figure EJ the spectrum is spread 
up to 1 ps, which corresponds to a large difference be- 
tween Ei and £3. Above that delay time, both averages 
coincide. At short times, an initial red shift is observed, 
whereas above 1 ps, a sudden blue shift occurs up to 1.3 
ps, followed by a final red shift: this panel, compared 
with figure m shows that the results are once again con- 
sistent with the expected time evolution of the plasmon. 
Even if well defined values of Ei and £"3 are hard to as- 
sign in the 300-800 fs delay range, one can see that the 
average values oscillate with a period of a few hundred 
fs which shrinks somewhat with time going on. We ex- 
pect that the mean plasmon frequency depends on the 
ionic structure. Thus we check the time evolution of the 
transverse and the total mean square radii, respectively 
(R±_ 2 ) and (R 2 ), see middle and top panels of figure 0] 
Note that after 1.2 ps, the radii are calculated for the two 
separate fission fragments relative to each separate cen- 
ter of mass. Before scission time, (R 2 ) shows a breathing 
although partially masked by the global extension of the 
cluster along its symmetry axis. The transversal radius 
(Rx ) also oscillates but is not exactly in phase with 
(R 2 ). Both radii present periods that shrink with in- 
creasing time. We note a relatively good correspondence 
between the minima of u> p i and the maxima of the radii, 
and vice versa for the maxima of w p i and the minima of 
the radii. The oscillations of the plasmon frequency are 
not dominated by either the transverse ionic motion or 
the longitudinal one, but seem to be a mixture of both. 
After scission, the plasmon evolution is clearly related 
to the radial oscillations since the extrema of u> p \ and 
(Ri 2 ) perfectly coincide. Thus this "measurement" of 
the plasmon frequency through the resonant ionization 
of the cluster gives access to the time evolution of the 
ionic structure, and to some extent to the potential en- 
ergy surface along the fission path. 

In this paper we have presented a scenario for pump 
and probe (P&P) analysis of fission of metal clusters. 
The scheme exploits the marked response of the Mie sur- 
face plasmon to laser fields and the fact that the plasmon 
resonance spectrum changes dramatically along the fis- 
sion path. The response is quantified in a measurable way 
in terms of ionization induced by the laser pulse. With 
a single-frequency P&P setup, one can measure the scis- 
sion time (that is, the final break up into fragments) when 
properly tuning the frequency. Scanning the frequency of 
the probe pulses would allow to resolve more details of 
the ionic geometries along the fission path. Finally, we 
ought to mention that a proper orientation of the clus- 



ter is required during the time span of the measurement 
time, i.e. for a few ps. This can be achieved by preparing 
preliminary a sample of clusters of same orientation, for 
instance, with the laser burning technique [23 . |23| ; and 
by cooling the beam to a temperature of about 10 K, so 
that the orientation is conserved for at least 2 ps. 
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